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We have carried out a Schrodinger functional (SF) calculation for the SU(3) lattice gauge theory 
with two flavors of Wilson fermions in the sextet representation of the gauge group. We find 
that the discrete beta function, which governs the change in the running coupling under a discrete 
change of spatial scale, changes sign when the SF renormalized coupling is in the neighborhood of 
g 2 = 2.0. The simplest explanation is that the theory has an infrared-attractive fixed point, but 
more complicated possibilities are allowed by the data. While we compare rescalings by factors of 
2 and 4/3, we work at a single lattice spacing. 
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I. INTRODUCTION 



Gauge theories with groups larger than SU(3), or with 
light fermions in representations higher than the funda- 
mental, are a staple of theories that go beyond the Stan- 
dard Model [l| . Among the mechanisms proposed to con- 
nect these theories to the Standard Model at low energies 
are technicolor 0, 0] and tumbling Q , with many asso- 
ciated variants. Both of these depend largely on weak- 
coupling pictures for their dynamics: a perturbative (3 
function to take technicolor from weak to strong coupling 
as the energy scale drops, and a most-attractive-channel 
argument for scale separation and selection of the con- 
densed channel in tumbling. Nonperturbative tests of 
these pictures are long overdue. While some lattice stud- 
ies have been carried out on SU(3) gauge theories with 
Nf > 3 fundamental flavors [a, la, LD, |8|, |f| , scant attention 
has been paid to more general alternatives where a richer 
set of phenomena may be sought [ll|, [H, [l3[ . 

A number of ideas focus on the behavior of the gauge 
theory's (3 function and the possibility of an IR-attractive 
fixed point [TJ, [l5[ . The infrared limit of the massless 
theory is then scale-invariant and probably conformal, 
devoid of confinement and of chiral symmetry break- 
ing [l6|, [ljj ■ Alternatively, the near-appearance of a fixed 
point makes the (3 function hover near the axis with- 
out crossing it; this is the scenario of "walking" pH. fli^. 
wherein enormous scale ratios are generated before con- 
finement finally sets in at large distances. 

Let us describe the perturbative picture of the theory 
we have studied. This is the SU(3) gauge theory with 
Nf fermions in the symmetric two-index representation 
[20l |2lj . which for this group is the sextet. Consider 
varying Nf upwards from zero. At first, the one- and 
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are both negative: this is like ordinary QCD. When Nf 
passes HI ~ 1.22, the two- loop coefficient in the (3 func- 
tion becomes positive and hence the two-loop (3 function 
acquires a zero at positive coupling g = g* . For the 
Nf = 2 theory studied in this paper, this zero is at the 
very large coupling g 2 ~ 10.4. As Nf grows g* becomes 
weaker, lending perha ps more credibility to the pertur- 
bative prediction [23|, |24| | . It should be kept in mind, 
though, that when one-loop and two-loop effects com- 
pete, there are usually similar-size contributions coming 
from higher orders in perturbation theory. An infrared 
fixed point generally requires nonperturbative confirma- 
tion. 

We have begun a study of the SU(3) lattice theory 
with Nf = 2 Wilson fermions in the sextet representa- 
tion. Since the (3 function gives the most direct approach 
to all these scenarios, we chose it as our first object of 
study. We apply the Schrodinger functional (SF) method 
HE ES HE EI, HI H HI S3, wherein we impose a 
background gauge field and calculate directly the running 
coupling as the scale of the background field is changed. 
Since we change the scale by a discrete amount, we ob- 
tain the discrete beta function (DBF), analogous to the 
usual (3 function. 

We find that the DBF of the massless theory crosses 
zero at g 2 ~ 2.0, far short of the perturbative prediction. 
If the full renormalization-group flow of the theory can 
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indeed be summed up by this single coupling constant, 
then this is an IR-attractive fixed point, implying scale 
invariance in the IR physics of the strictly massless lat- 
tice theory defined in the fixed point's catchment basin. 
At the same time, the zero of the DBF allows for more 
complex possibilities stemming from more complex RG 
flows, as discussed below. 

In order to judge the significance of the zero of the 
DBF, as well as to study the physics of the theory in its 
vicinity, we have also calculated observables connected 
with the qq potential and with chiral symmetry. While 
our presentation here is brief, we wish to place it in its 
proper context. A simple and general argument shows 
that the formation of any bound states made out of light 
fermions in a gauge theory (without scalars) necessitates 
the spontaneous breaking of chiral symmetry 33]. Thus, 
if the infrared limit is a confining theory, chiral symme- 
try must be broken spontaneously at an en ergy scale at 
least as high as the confinement scale [33j, [34J. 1 If the 
two scales are indeed separated, the intermediate region 
breaks chiral symmetry but does not show confinement; 
and an entirely nonconfining theory can still break chi- 
ral symmetry spontaneously. A conformal theory, on the 
other hand, is inconsistent with spontaneously broken 
chiral symmetry since the breaking creates a mass scale. 



II. SCHRODINGER FUNCTIONAL METHOD 

In order to introduce the SF running coupling g{L), 
we begin with the gauge theory defined in a small Eu- 
clidean box of volume L 4 . We can consistently choose the 
coupling in this volume to be small; asymptotic freedom 
ensures that there is only one effective coupling, that it 
runs with the perturbative /3 function, and hence that 
at yet smaller distances, the coupling is even smaller. 
In determining the running coupling non-perturbatively, 
virtually any observable can in principle be used to ex- 
tract it. For consistency, we require that a perturbative 
calculation of the same observable will indeed reproduce 
the running coupling at this small scale. 

The SF definition of the running coupling g(L) is an 
application of the background field method. Consider a 
background field calculation in the classical field strength 
Fftv/g. If by construction the only distance scale that 
characterizes the background field is L, the n-loop effec- 
tive action T = — log Z gives the running coupling via 



T = g{L)- 2 S 



where 
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1 Casimir scaling is the most popular mechanism to explain how 
chiral symmetry could be broken at an energy scale much higher 
than the confinement scale fiol . Illl . Il2t |35| . |36| . 



g(L) is the result of integrating the n-loop f3 function. 

The SF defines the background field by imposing 
Dirichlet boundary conditions at t = and t = L. The 
boundary values are chosen such that the classical action 
has a unique, nontrivial minimum, and the configuration 
at this minimum is F^jg 25]. The effective action Y 
is then calculated and compared via Eq. ^ to Sy M , 
which in turn is to be evaluated for the classical field 
that minimizes it with the given boundary conditions. If 
the effective action is calculated nonperturbatively, this 
procedure gives a nonperturbative definition of g{L). 

At short distances, the static potential between color 
sources is a Coulomb potential. The scale dependence of 
g 2 provides a small correction. We postulate that this is 
true when the volume is small enough that confinement 
is not evident. In large volumes, on the other hand, the 
static potential can be qualitatively different; the notion 
of a unique effective coupling that depends only on the 
overall scale may no longer be tenable. The upshot is 
that one must be cautious in drawing conclusions from 
the running of a single coupling constant. We bear this 
in mind as we take Eq. (UJ) to define of g(L) beyond per- 
turbation theory. 

The continuum framework carries over to the lattice 
with minor adaptations. A technical obstacle is that 
Monte Carlo methods do not allow for the direct com- 
putation of the effective action. This is solved by consid- 
ering a family of gauge-field boundary values that depend 
on a continuous parameter r/. By differentiating Eq. (J2|) 
we obtain 
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The derivative of T gives an observable quantity, while 
K is just a number [261 ]. 

Our lattice theory is defined by the single-plaquette 
gauge action and a Wilson fermion action with added 
clover term [13] ■ Our choice of Nf = 2 ensures positiv- 
ity of the fermion determinant, which in turn allows use 
of the standard hybrid Monte Carlo algorithm. We es- 
chew perturbative corrections [291 ] to the parameters and 
operators in this exploratory study, but we do modify 
the clover term's coefficient via self-consistent tadpole 
improvement, csw = l/ u o- The quantum correction 
csw — 1 is known to be nonperturbative and large for 
fundamental fermions (3lj . where it serves to bring the 
theory closer to the continuum limit; the tadpole esti- 
mate supplies a large portion of the correction. 

The SF boundary conditions take the form of fixed, 
spatially constant values for the spacelike links Ui on the 
top and bottom layers of the lattice. These links enter 
into the gauge plaquette and into the clover term of the 
fermion action. Thus the r\ derivative of the effective 
action is given by 
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where Dp is the complete Wilson-clover fermion action. 
(The appearance of D F Dp indicates that Nf = 2; we 
evaluate the functional trace with a noisy estimator 1321. ) 
The boundary fields are chosen as described in Ref. [3l| ; 
the parameter rj enters linearly into phase angles so the 
derivatives on the right-hand side of Eq. (|7|) are imple- 
mented by putting the appropriate fixed values in place 
of the boundary links [26j. With these boundary values 
the coefficient 2 K = 37.7. We also impose twisted spatial 
boundary conditions on the fermion fields as suggested in 
Ref. [H, iMx + L) = exp(iff)ip(x), with 6 = tt/5 on all 
three axes |32j |. 

Equation ^ defines the running coupling at any given 
length scale L, which we take to be the linear size of the 
lattice. The discrete beta function gives the change in 
K I g 2 when L is multiplied by n. Defining u = K/g 2 (L), 
we write 

B(u,n) = ^—- U , (8) 
g\nL) 

which is the counterpart of Eq. Jl}. 3 We use a lat- 
tice approximation of Eq. ([8]), which introduces an im- 
plicit dependence on the ultraviolet cutoff: We calcu- 
late B(u,n — 2) with lattice spacing a = L/4, and 
B(u 7 n = 4/3) with a = L/6. 

ILL LATTICE CALCULATION OF THE 
DISCRETE BETA FUNCTION 

We study herein the scaling of only the massless theory. 
In the case of Wilson fermions, where the quark mass is 
unprotected against additive renormalization, this means 
fixing the hopping parameter to its critical value, k = n c , 
at each value of the bare lattice coupling (3 = 6/g^. It 
is easy to locate k c (/3) when SF boundary conditions are 
used, since these boundary conditions (and the spatial 
twists) limit the condition number of the fermion matrix 
even when there is no mass, and hence we can carry out 
simulations precisely at k c .A straightforward way of find- 
ing k c is to calculate the quark mass as defined by the 
lattice approximation of an axial Ward identity (AWI), 

_ 1 d 4 (A b (t)O b (0)) 
mq ~ 2 {P b (t)O b (0)) ■ ( ' 

Here A\(t) — 'tp^j^tp is the time component of the 
local axial vector current with flavor 6, taken at zero 
spatial momentum on the time slice t; P b {t) is the lo- 
cal pseudoscalar density. The operator O b (0) is defined 
by introducing spatially constant (Grassmann) Dirich- 
let boundary conditions for the fermions, differentiating 



2 This is the value for a lattice with 8 4 sites. There is a tiny 
dependence on the lattice size, which we safely neglect. 

3 In the SF literature it is customary to define the step scaling 
function cr(v) = g 2 (2L), with v = g 2 (L). 



TABLE I: k c and uo for L — 4a with SF boundary conditions. 
Linear interpolation may be used safely between (3 — 5.0 
and 5.5 and between j3 = 5.5 and 6.0. 






Kc 
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5.0 


.1723 


.875 


5.5 


.1654 


.887 


6.0 


.1610 


.900 


7.0 


.1536 


.916 


8.0 


.1486 


.928 



with respect to these boundary values, and finally setting 
them to zero [29[ . The derivative in Eq. © is a symmet- 
ric difference evaluated about t = L/2, the center of the 
lattice. Equation ([9]) neglects multiplicative renormaliza- 
tion of the currents, but this is unimportant since we only 
use it to locate k c by demanding m q = 0. In addition, 
we neglect the mixing of A b ^ with d^P b ', which is known 
to be a small effect for fundamental fermions [2!| [3l[ . 

We list the values of k c and the tadpole coefficient i*o in 
Tabled] While these values were determined for L — 4a, 
we find a small movement in m q (0.01-0.04) when going 
to L = 6a and then L — 8a. Calculation of the DBF de- 
mands keeping the ultraviolet cutoff fixed as the volume 
is changed, and thus /?, K, and uq must remain unchanged 
when comparing different lattices. The undesirable shift 
in uiq is a discretization error; it can only be reduced by 
going to larger pairs of lattices. 

Beginning with the scale factor n = 2, we thus cal- 
culate the running coupling K/ g 2 for given bare cou- 
pling P at the critical hopping parameter k — k c (J3), 
first on a lattice with 4 4 sites, which defines the scale 
L = 4a, and then at the same ((3, k) on an 8 4 lat- 
tice, which gives the scale 2L. We show the results in 
Fig. [TJ At large (3 (corresponding to the perturbative 
regime, large u) the DBF agrees with the one-loop re- 
sult B(u,2) = -[X6i/(167r 2 )]21og2 ~ -1.43; the dashed 
curve shows the two-loop result. It is plain that the DBF 
departs from two-loop perturbation theory and crosses 
the axis in the neighborhood of K/g 2 — 19, or g 2 ~ 2.0. 
The corresponding bare coupling is (3 ~ 5.6. 

For an indication of the dependence of the result on 
the lattice spacing a = L/4, we have carried out a par- 
allel calculation of B(u,A/3) with a = L/6. This means 
calculating K/ g 2 on a lattice with 6 4 sites, at the same 
(3 values as above, 4 and comparing to the results on 8 4 
sites. The first lattice defines the scale L while the sec- 
ond now gives the scale 4L/3. This is economical because 
the same 8 4 data are used in this calculation as in the 
preceding. 

These results are shown in Fig. [2] along with the two- 
loop prediction. The latter is smaller than that for Fig.[T] 
because the scale factor is 4/3 rather than 2. In the 



4 We use the same values of re = k c and uq as above. 
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u = K/g 2 (4 4 ) 



FIG. 1: Discrete beta function for the scale transformation 
L —* 2L, as defined in Eq. (JSJ. The lattice spacing is fixed 
such that L — 4a. The dashed curve is the two- loop result. 
The data points are calculated at bare couplings (left to right) 
(3 = 5.2 to 6.0 by 0.1, and then (3 = 7.0 and 8.0. Horizontal 
error bars are the size of the plotted symbols. 



lattice data, the crossing of zero (again near g 2 = 2.0) 
is evident. This comparison of course does not take the 
place of a real scaling study. 



IV. NONCONFINEMENT AND CHIRAL 
SYMMETRY 

We have also carried out calculations on lattices of 
volume 8 3 x N t for N t = 8 and 12 with ordinary 
(anti-)periodic boundary conditions, at several values of 
k < k c in the same range of (3. These calculations show 
that (3 = 5.6, with k = k c , lies well on the weak-coupling 
side of the finite-volume crossover associated with con- 
finement physics: The average Polyakov loop is large and 
the static qq potential is entirely Coulombic in the range 
allowed by the volume. This means that the lattices of 
size L = 8 are not large enough to contain the physics of 
confinement. As discussed above, this is essential for in- 
terpreting g{L) as an effective coupling that characterizes 
the theory. 

By QCD standards, a coupling g 2 /Ait ~ 0.16 is per- 
turbative. This is in line with the lack of confine- 
ment in our 8 4 lattices. Nevertheless, this coupling may 
not be weak when considering the possibility of sponta- 
neous chiral symmetry breaking for higher-representation 
fermions. Indeed, a quenched study with color-sextet 
fermions found that the finite-temperature chiral transi- 
tion occurs at f3 ~ 7.8 with Nt = 4, a far larger (3 value 
than that of the deconfinemcnt transition 10]. With our 
dynamical fermions, however, things are different. On 
lattices with 8 2 x 12 x N t sites, with N t — 8, we find the 
chiral restoration crossover at k close to, but definitely 
below K c ((3). Thus the coupling at the zero of the DBF 
doesn't break chiral symmetry either. 



DISCUSSION 
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FIG. 2: Discrete beta function for the scale transformation 
L — * 4L/3. The lattice spacing is fixed such that L — 6a. 
The dashed curve is the two-loop result. The data points are 
calculated at the same bare couplings as in Fig. [1] Horizontal 
error bars are the size of the plotted symbols. 



The simplest explanation of the zero of the DBF is 
that the strictly massless theory with two color-sextet 
fermions has an infrared-attractive fixed point, leading to 
conformal physics without chiral symmetry breaking and 
without confinement. Accepting this conclusion, a major 
goal of further lattice studies would be to understand the 
behavior of the same theory at nonzero fermion mass. In 
a truly conformal theory, the introduction of a fermion 
mass provides the only scale; that scale can sometimes 
play the role of an effective ultraviolet cutoff, and some- 
times the role of an infrared cutoff. This is an unfamiliar 
territory that contains many interesting questions. 

Our results clearly do not allow for a continuum ex- 
trapolation. More lattice volumes are needed, and it is 
also desirable to compute the DBF for more scale ratios. 
Such detailed information would allow the reconstruction 
of a continuous (3 function. Appelquist, Fleming, and 
Neil [|[ have recently presented an extensive SF analy- 
sis of the SU(3) gauge theory with Nf — 8 and 12 fla- 
vors of fundamental fermions. Their results are based on 
larger statistics as well as a detailed continuum extrapo- 
lation. This allows them to conclude with certainty that 
the Nf — 12 theory has an infrared fixed point. 

Our results for the DBF of the color-sextet theory 
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do not preclude more elaborate scenarios whereby the 
dynamics generates new effective degrees of freedom at 
some nonperturbative scale. This would typically lead 
to new relevant and/or marginal couplings; our single- 
parameter DBF would result from projecting the multi- 
parameter RG flow into a one-dimensional subspace. 

A concrete scenario in this direction was proposed by 
one of us in Ref. [38|. One supposes that chiral sym- 
metry breaks spontaneously when the interaction is still 
Coulombic. This allows for the existence of colored ex- 
citations. The fermion condensate is then polarized in 
response to an applied color field. As we move down 
in energy scale across the chiral transition, the effective 
gauge coupling is screened by the collective response of 
the fermion condensate. This effect, if found, would be 
the relativistic counterpart of the familiar dielectric po- 
larization. The scenario is described by an effective La- 
grangian that indeed contains new couplings 38]. For a 
short while, the running of the gauge coupling reverses 
its direction from that of an asymptotically free theory; 
a discrete sampling of this running can result in a zero of 
the DBF, imitating a fixed point. 

Whatever the scenario, our DBF flatly contradicts con- 
tinuum perturbative estimates. If confirmed by a true 
scaling study, it would show that this theory is not 



like QCD with a small number of fundamental fermions, 
where the beta function is always negative; and it is not 
like theories where the fermions condense and decouple 
from the gauge fields, since in that case the beta function 
would behave as if Nf = 0, becoming even more negative 
than that of the fully coupled theory. The simplest ex- 
planation of our DBF is still an infrared-attractive fixed 
point, at an unexpectedly weak coupling. This would 
place the massless theory squarely in the conformal win- 
dow: It would prevent confinement and the spontaneous 
breaking of chiral symmetry, and remove the theory from 
the list of candidates for walking. 

We will present a detailed study of the physics of this 
model, so different from QCD, in a forthcoming paper. 
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